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ABSTRACT 


We propose the quantum simulation of the quantum Rabi model in all parameter regimes by means of detuned bichromatic 
sideband excitations of a single trapped ion. We show that current setups can reproduce, in particular, the ultrastrong and deep 
strong coupling regimes of such a paradigmatic light-matter interaction. Furthermore, associated with these extreme dipolar 
regimes, we study the controlled generation and detection of their entangled ground states by means of adiabatic methods. Ion 
traps have arguably performed the first quantum simulation of the Jaynes-Cummings model, a restricted regime of the quantum 
Rabi model where the rotating-wave approximation holds. We show that one can go beyond and experimentally investigate the 
quantum simulation of coupling regimes of the quantum Rabi model that are difficult to achieve with natural dipolar interactions. 


Introduction 

The quantum Rabi model (QRM) describes the most fundamental light-matter interaction involving quantized light and 
quantized matter. It is different from the Rabi model/ where light is treated classically. In general, the QRM is used to describe 
the dipolar coupling between a two-level system and a bosonic field mode. Although it plays a central role in the dynamics of 
a collection of quantum optics and condensed matter systems,^ such as cavity quantum electrodynamics (CQED), quantum 
dots, trapped ions, or circuit QED (cQED), an analytical solution of the QRM in all coupling regimes has only recently been 
proposed.^’ ^ In any case, standard experiments naturally happen in the realm of the Jaynes-Cummings (JC) model,^ a solvable 
system where the rotating-wave approximation (RWA) is applied to the QRM. Typically, the RWA is valid when the ratio 
between the coupling strength and the mode frequency is small. In this sense, the JC model is able to correctly describe most 
observed effects where an effective two-level system couples to a bosonic mode, be it in more natural systems as CQED,^“^ or 
in simulated versions as trapped ions^’ and cQED.^^’ However, when the interaction grows in strength until the ultrastrong 
coupling (USC)^^“^^ and deep strong coupling (DSC)^^’ regimes, the RWA is no longer valid. While the USC regime happens 
when the coupling strength is some tenths of the mode frequency, the DSC regime requires this ratio to be larger than unity. In 
such extreme cases, the intriguing predictions of the full-fiedged QRM emerge with less intuitive results. 

Recently, several systems have been able to experimentally reach the USC regime of the QRM, although always closer 
to conditions where perturbative methods can be applied or dissipation has to be added. Accordingly, we can mention the 
case of circuit QED,^"^’^^ semiconductor systems coupled to metallic microcavities,split-ring resonators connected to 
cyclotron transitions,^^ or magnetoplasmons coupled to photons in coplanar waveguides.The advent of these impressive 
experimental results contrasts with the difficulty to reproduce the nonperturbative USC regime, or even approach the DSC 
regime with its radically different physical predictions.^^’ Nevertheless, these first achievements, together with recent 

theoretical advances, have put the QRM back in the scientific spotlight. At the same time, while we struggle to reproduce the 
subtle aspects of the USC and DSC regimes, quantum simulators'^ may become a useful tool for the exploration of the QRM 
and related models. 

Trapped ions are considered as one of the prominent platforms for building quantum simulators.^ ^ In fact, the realization 
and thorough study of the JC model in ion traps, a model originally associated with CQED, is considered a cornerstone in 
physics.This is done by applying a red-sideband interaction with laser fields to a single and may be arguably 

presented as the first quantum simulation ever implemented. In this sense, the quantum simulation of all coupling regimes 
of the QRM in trapped ions would be a historically meaningful step forward in the study of dipolar light-matter interactions. 
In this article, we propose a method that allows the access to the full-fledged QRM with trapped-ion technologies by means 


of a suitable interaction picture associated with inhomogeneously detuned red and blue sideband excitations. Note that, in 
the last years, bichromatic laser fields have been successfully used for different purposes.In addition, we propose an 
adiabatic protocol to generate the highly-correlated ground states of the USC and DSC regimes, paving the way for a full 
quantum simulation of the experimentally elusive QRM. 


Results 

Single atomic ions can be confined using radio-frequency Paul traps and their motional quantum state cooled down to its ground 
state by means of sideband cooling techniques. In this respect, two internal metastable electronic levels of the ion can play 
the role of a quantum bit (qubit). Driving a monochromatic laser field in the resolved-sideband limit allows for the coupling of 
the internal qubit and the motional mode, whose associated Hamiltonian reads 

H = va^a-\-Q.{G^ + G~) (^exp{i[ri{a-{-a^) — coit + 0/]} -f exp{—/[r 7((2 + a^) — coit + 0/]}^. 


Here, and a are the creation and annihilation operators of the motional mode, cr+ and G are the raising and lowering Pauli 
operators, v is the trap frequency, COq is the qubit transition frequency, Q. is the Rabi coupling strength, and 7] = is 

the Lamb-Dicke parameter, where k is the component of the wave vector of the laser on the direction of the ions motion and 
m the mass of the ion;^ (0/ and 0/ are the corresponding frequency and phase of the laser field. For the case of a bichromatic 
laser driving, changing to an interaction picture with respect to the uncoupled Hamiltonian, F/q = + va^a and applying an 

optical RWA, the dynamics of a single ion reads^ 







( 1 ) 


with a{t) = ae and We will consider the case where both fields are off-resonant, first red-sideband (r) and 

first blue-sideband (b) excitations, with detunings 5r and 5^, respectively. 


(Or = (OQ-V^-5r, U)^ = (00 + V + 4. 


In such a scenario, one may neglect fast oscillating terms in Eq. (1) with two different vibrational RWAs. We will restrict 
ourselves to the Lamb-Dicke regime, that is, we require that a) 1. This allows us to select terms that oscillate 

with minimum frequency, assuming that weak drivings do not excite higher-order sidebands, 5n^kln ^ V for n = r^b. These 
approximations lead to the simplified time-dependent Hamiltonian 

+H.c., (2) 

where we consider equal coupling strengths for both sidebands, D = = D/,. Equation 2 corresponds to the interaction picture 

Hamiltonian of the Rabi Hamiltonian with respect to the uncoupled Hamiltonian F/q = \ ^ 8r)Gz+ \ — 8r)a^a, 

co^ 

Hqrm = -^(y^ + (0^a^a + ig{(y^ -(y~){a + a^), (3) 

with the effective model parameters given hy 

4 = - ^(5. + 4), = ^(5.- 4), g = ^, (4) 

where the qubit and mode frequencies are represented by the sum and difference of both detunings, respectively. The tunability 
of these parameters permits the study of all coupling regimes of the QRM via the suitable choice of the ratio g/(0^. It is 
important to note that all realized interaction-picture transformations, so far, are of the form aa^a + jScj^. This expression 
commutes with the observables of interest, {cr^, \n){n\,a'^a}, warranting that their experimental measurement will not be 
affected by the transformations. 


Accessible regimes 

The quantum Rabi model in Eq. (3) will show distinct dynamics for different regimes, which are defined by the relation among 
the three Hamiltonian parameters: the mode frequency (0^, the qubit frequency U)q , and the coupling strength g. 

We first explore the regimes that arise when the coupling strength is much weaker than the mode frequency \co^\- 
Under such a condition, if the qubit is close to resonance, jco^l ^ |(0 q |, and | + (Oq | | (0^ — (Oq | holds, the RWA can be 


2/9 





applied. This implies neglecting terms that in the interaction picture rotate at frequency C0^-\-C0q, leading to the JC model. This 
is represented in Fig. 1 by the region 1 in the diagonal. Notice that these conditions are only possible if both the qubit and the 
mode frequency have the same sign. However, in a quantum simulation one can go beyond conventional regimes and even reach 
unphysical situations, as when the qubit and the mode have frequencies of opposite sign. In this case, \(0^ — COq\ ^ | + (Oq | 

holds, see region 2, and we will be allowed to neglect terms that rotate at frequencies \co^ — COq\. This possibility will give rise 

to the anti-Jaynes Cummings (AJC) Hamiltonian, co^a^a + ig{a^a^ — a~a). It is noteworthy to mention that, 

although both JC and AJC dynamics can be directly simulated with a single tuned red or blue sideband interaction, respectively, 
the approach taken here is fundamentally different. Indeed, we are simulating the QRM in a regime that corresponds to such 
dynamics, instead of directly implementing the effective model, namely the JC or AJC model. 

If we depart from the resonance condition and have all terms rotating at high frequencies {| |, | (Oq |, | + (Oq |, | — 

^0 1} ^ region 3, for any combination of frequency signs, the system experiences dispersive interactions governed by a 
second-order effective Hamiltonian. In the interaction picture, this Hamiltonian reads 


Heff = g^ 


\e){e\ 

co^ - co^ 


\ 8 ){ 8 \ 

CO^-^COq 


2co^ 


{co^ ^ co^){co^ - co^) 


a^ao. 


(5) 


inducing AC-Stark shifts of the qubit energy levels conditioned to the number of excitations in the bosonic mode. 

The use regime is defined as 0.1 ^g/co^ ^ 1, with perturbative and nonperturbative intervals, and is represented in Fig. 1 
by region 4. In this regime, the RWA does not hold any more, even if the qubit is in resonance with the mode. In this case, 
the description of the dynamics has to be given in terms of the full quantum Rabi Hamiltonian. For > 1, we enter into 
the DSC regime, see region 5 in Fig. 1, where the dynamics can be explained in terms of phonon number wave packets that 
propagate back and forth along well defined parity chains. 

In the limit where co^ = 0, represented by a vertical centered red line in Fig. 1, the quantum dynamics is given by the 
relativistic Dirac Hamiltonian in 1-^1 dimensions. 


Hd = mc^Gz + cpGx, ( 6 ) 

which has been successfully implemented in trapped ions,^^’^^ as well as in other platforms.^^’"^^ 

Moreover, an interesting regime appears when the qubit is completely out of resonance and the coupling strength is small 
when compared to the mode frequency, COq ^ 0 and g <C | |. In this case, the system undergoes a particular dispersive 

dynamics, where the effective Hamiltonian becomes a constant. Consequently, the system does not evolve in this region that we 
name as decoupling regime, see region 6 in Fig. 1. The remaining regimes correspond to region 7 in Fig. 1, associated with the 
parameter condition \ cOq\ ^ g <C | |. 

The access to different regimes is limited by the maximal detunings allowed for the driving fields, which are given by the 
condition 5^^^ V, ensuring that higher-order sidebands are not excited. The simulations of the JC and AJC regimes, which 
demand detunings \5r^t>\ < |m^| + I^q |, are the ones that may threaten such a condition. We have numerically simulated the 
full Hamiltonian in Eq. (1) with typical ion-trap parameters: V = 271 x 3MHz, Q = 27r x 68kHz and rj = 0.06,^^ while the 
laser detunings were 5^ = —271 x 102kHz and 5r = 0, corresponding to a simulation of the JC regime with g/(o^ = 0.01. The 
numerical simulations show that second-order sideband transitions are not excited and that the state evolution follows the 
analytical JC solution with a fidelity larger than 99% for several Rabi oscillations. This confirms that the quantum simulation 
of these regimes is also accessible in the lab. We should also pay attention to the Lamb-Dicke condition rj ^/{a^ <C 1, as 
evolutions with an increasing number of phonons may jeopardize it. However, typical values like T] = 0.06 may admit up to 
some tens of phonons, allowing for an accurate simulation of the QRM in all considered regimes. 

Regarding coherence times, the characteristic timescale of the simulation will be given by ^char ~ simulator, 

g = ^, such that ^char = typical values of r] = 0.06 — 0.25 and of ^l/2;r = 0 — 500 kHz, the dynamical timescale of 

the system is of milliseconds, well below coherence times of qubits and motional degrees of freedom in trapped-ion setups. 


State preparation 

The ground state |G) of the QRM in the JC regime (g <C co^) is given by the state |g,0), that is, the qubit ground state, |g), and 
the vacuum of the bosonic mode, |0). It is known that |g, 0) will not be the ground state of the interacting system for larger 
coupling regimes, where the contribution of the counter-rotating terms becomes important.As seen in Fig. 2, the ground state 
of the USC/DSC Hamiltonian is far from trivial,^ essentially because it contains qubit and mode excitations, {G\a'^a\G) > 0. 

Hence, preparing the qubit-mode system in its actual ground state is a rather difficult state-engineering task in most 
parameter regimes, except for the JC limit. The non analytically computable ground state of the QRM has never been observed 
in a physical system, and its generation would be of theoretical and experimental interest. We propose here to generate the 
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ground state of the USC/DSC regimes of the QRM via adiabatic evolution. Figure 3 shows the fidelity of the state prepared 
following a linear law of variation for the coupling strength at different evolution rates. When our system is initialized in the 
JC region, achieved with detunings 5^ = 0 and 15/, | ^ g, it is described by a JC Hamiltonian with the ground state given by 
|G) = |g, 0). Notice that the g/d)^ ratio can be slowly turned up, taking the system adiabatically through a straight line in the 
configuration space to regions 4-5."^^ This can be done either increasing the value of g by raising the intensity of the driving, or 
decreasing the value of co^ by reducing the detuning 15^ |. The adiabatic theorem"^^ ensures that if this process is slow enough, 
transitions to excited states will not occur and the system will remain in its ground state. As expected, lower rates ensure a 
better fidelity. 

Once the GS of the QRM is generated, one can extract the populations of the different states of the characteristic parity 
basis shown in Fig. (2). To extract the population of a specific Fock state, one would first generate a phonon-number dependent 
ac-Stark shift."^^ A simultaneous transition to another electronic state will now have a frequency depending on the motional 
quantum number. By matching the frequency associated with Fock state n, we will flip the qubit with a probability proportional 
to the population of that specific Fock state. This will allow us to estimate such a population without the necessity of 
reconstructing the whole wave function. 

Discussion 

We have proposed a method for the quantum simulation of the QRM in ion traps. Its main advantage consists in the accessibility 
to the USC/DSC regimes and the convenient switchability to realize full tomography, outperforming other systems where the 
QRM should appear more naturally, such as cQED."^^’^^ In addition, we have shown how to prepare the qubit-mode system in 
its entangled ground state through adiabatic evolution from the known JC limit into the USC/DSC regimes. This would allow 
for the complete reconstruction of the QRM ground state, never realized before, in a highly controllable quantum platform as 
trapped ions. The present ideas are straightforwardly generalizable to many ions, opening the possibility of going from the 
more natural Tavis-Cummings model to the Dicke model. In our opinion, the experimental study of the QRM in trapped ions 
will represent a significant advance in the long history of dipolar light-matter interactions. 
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Figure 1. Configuration space of the QRM. (1) JC regime: g ^ IcOq |} and \co^ — C0 q\ + cOq |. (2) AJC regime: 

g <C \(0q\} and |(0^ — cOq | > \co^ + C0q\. (3) Two-fold dispersive regime: g < {|tiL)^|, IcOq |, — cOq |, \co^ + ^o\}- (4) 

use regime: Im^l < lOg, (5) DSC regime: Im^l < g, (6) Decoupling regime: \cOq\ g 1(5)^|. (7) This intermediate regime 
(\cOq\ ^ g |ciL)^|) is still open to study. The (red) central vertical line corresponds to the Dirac equation regime. Colours 
delimit the different regimes of the QRM, colour degradation indicates transition zones between different regions. All the areas 
with the same colour correspond to the same region. 


7/9 








0.2 


0.1 


0 el g2 e3 ^4 e5 ^6 el gS e9 glO ellgl2 el3gU elb 

Figure 2. State population of the QRM ground state. We plot the case of g/co^ = 2, parity /? = +!, and corresponding parity 
basis {|g,0), |e, 1), |g,2), |e, 3),...}. Here, p is the expectation value of the parity operator P = and only states 

with even number of excitations are populated. We consider a resonant red-sideband excitation (5^ = 0), a dispersive 
blue-sideband excitation (5^/2;r = —11.31kHz), and g = — 5^, leading to the values co^ = (Oq = —5^/2 and gjco^ = 2. 
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Figure 3. Fidelity of the adiabatic evolution. Let us assume that the system is initially prepared in the JC ground state |g,0), 
that is, when g co^. Then, the coupling is linearly chirped during an interval At up to a final value gf, i. e., g{t) = gftjAt. 
For slow changes of the laser intensity, the ground state is adiabatically followed, whereas for non-adiabatic processes, the 
ground state is abandoned. The instantaneous ground state \G{t)) is computed by diagonalizing the full Hamiltonian at each 
time step, while the real state of the system \ \i/{t)) is calculated by numerically integrating the time dependent Schrodinger 
equation for a time-varying coupling strength g{t). For the simulation, a ^^Ca+ ion has been considered with parameters: 

V = 2;r X 3MHz, 5, = 0, 4 = -6 x IQ-^y, rj = 0.06 and % = 2;r x bSkHz.^i 
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